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Abstract 

The current-voltage characteristics in Langmuir-Blodgett monolayers of 7-hexadecylquinolinium 
tricyanoquinodimethanide (C16H33Q-3CNQ) sandwiched between Al or Au electrodes is calculated, 
combining ab initio and self-consistent tight binding techniques. The rectification current depends 
on the position of the LUMO and HOMO relative to the Fermi levels of the electrodes as in the 
Aviram-Ratner mechanism, but also on the profile of the electrostatic potential which is extremely 
sensitive to where the electroactive part of the molecule lies in the monolayer. This second effect 
can produce rectification in the direction opposite to the Aviram-Ratner prediction. 
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I. INTRODUCTION 



Molecular electronics has gained large-scale interest in recent years. A promising sys- 
tem is the electrical rectifier based on suitably engineered molecules. At the origin of this 
idea, Aviram and Ratner- proposed in 1974 to use D-cr-A molecules, where D and A are, 
respectively, an electron donor and an electron acceptor, and a is a covalent "sigma" bridge. 
Electrical rectification was observed only recently, in particular in Langmuir-Blodgett (LB) 
multilayers or monolayers of 7-hexadecylquinolinium tricyanoquinodimethanide (C16H33Q- 
3CNQ, 1, Fig. [1]) sandwiched between metallic electrodes^-. Even if these results represent 
an important progress to achieve molecular electronics, the physical mechanism responsible 
for the rectification is not clear. One critical issue is to know if the Aviram-Ratner model can 
be applied to 1, because it is a D(+)-7r-A(-) molecule-. On the theoretical side, these molec- 
ular diodes are complex systems, characterized by large and inhomogeneous electric fields, 
which result from the molecular dipoles 5 in the layer, the applied bias, and the screening 
induced by the molecules themselves and the metallic electrodes. A theoretical treatment of 
these effects is presently lacking, as it requires a self-consistent resolution of the quantum- 
mechanical problem, including the effect of the applied bias on the electronic structure. Our 
aim in this paper is to present such a theory, and its application to the systems experimen- 
tally studied in Refs.-^. Combining ab initio and semi-empirical calculations, we show that 
the direction of easy current flow (rectification current) depends not only on the position 
of the HOMO and LUMO, relative to the Fermi levels of the metal electrodes before bias 
is applied, but also on their shift after the bias is applied: this situation is more complex 
than the Aviram-Ratner mechanism, and can provide a rectification current in the opposite 
direction. We calculate that the electrical rectification results from the asymmetric profile 
of the electrostatic potential across the system. We obtain that the Ci6H 33 tail plays an 
important role in this asymmetry, and we predict a more symmetric I{V) curve in the case 
of molecules with a small alkyl chain. Quite generally, our work emphasizes the importance 
of the electrostatic potential profile in a molecular system and suggests that this profile 
could be chemically designed to build new devices. 
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II. ELECTRONIC STRUCTURE AND POTENTIAL 

To calculate the electronic structure and the geometry of the free molecules (gas phase), 
we have performed ab initio calculations, using the DMOL code 7 in the local density approx- 
imation (LDA) with the functional of Ref.-, and in the generalized gradient approximation 
(GGA) of Ref. 9 . The two approximations give very close results. We use a double numeri- 
cal basis set (two atomic orbitals for each occupied orbital in the free atom), together with 
polarization functions (2p for H, and 3d for N and C). As the molecular diode is a complex 
system, which cannot be treated in LDA, we have developed self-consistent tight binding 
(TB) calculations, as described in Ref.—. We start from the GGA optimized geometries. C 
and N atoms are represented by one s and three p atomic orbitals, and H atoms by one s 
orbital. The diagonal terms of the Hamiltonian matrix are charge-dependent^ 1 ^ (in atomic 
units): 

Hia )ia = H® a ia - i ^ J = - eV ext (i) (1) 

where a is the orbital index, i and j are atomic indices (j runs over all the atoms of 
the system), Qj is the net charge on atom j, Rij is the distance between atoms % and j, 
Uq is the intra-atomic Coulomb energy and the quantities Hf a ia define the s and p atomic 
levels. V ext (i) is the potential induced by the charges on the electrodes. It is calculated 
using a self-consistent solution of the Schrodinger and Poisson equations, with the boundary 
condition imposed by the applied bias. The second term in eqn. [T] is evaluated using the 
method of Ewald^ 2 -, taking into account the two-dimensional periodicity of the LB film. The 
molecules are ordered on a 9Ax 6.5A rectangular lattice, corresponding to a molecular area 
of 58.5 A 2 , close to the experimental value of ~ 50 A 2 -. The molecules are inclined by 45°, 
as measured experimentally^. Chemical interactions between molecules are supposed to be 
negligible. The LB film is sandwiched between two planar Al or Au electrodes, which we put 
at 3.2 A from the edges of the molecule, a reasonable value for molecules in weak interaction 
with a metallic surface^. We applied our calculations to 1, but 3 also to molecules with a 
smaller alkyl chain or without alkyl chain (Q-3CNQ, 2) for reasons which will become clear 
later. As expected, the electronic structure cose to the gap does not depend on the length 
of the alkyl chain for the free molecules. We predict a large dipole moment of 21 D in LDA 
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and 27 D in TB. The general results in LDA and in TB are quite similar, and agree with 
recent calculations^. Fig. [I] describes the original Aviram-Ratner mechanism^ for a D-a-A 
molecule. The electronic states are supposed to be totally localized either on the D side or 
on the A side. The HOMO(D) and LUMO(D) are high in energy, compared to, respectively, 
the HOMO (A) and LUMO(A). Therefore, a current can be established at relatively small 
positive bias, such that the Fermi level at the A side is higher than LUMO(A), and the 
Fermi level at the D side is lower than HOMO(D), providing that the electrons can tunnel 
inelastically through the a bridge. Thus we expect an asymmetric I(V) curve with an 
important rise of the current at positive bias above a given threshold (the electrode at A 
side is grounded, as in the experiments^^). In the case of our D(+)-7r-A(-) molecules 
(1 and 2), LDA and TB calculations show that the HOMO and the LUMO are strongly 
delocalized on the whole molecule. Even if the HOMO is more localized on the A part and 
the LUMO on the D part, the extension of the wave functions on the opposite side exceeds 
30%. Thus the two parts of the molecules cannot be decoupled, as in the original Aviram- 
Ratner model. This is obviously due to the 7r bridge between the donor and acceptor sites. 
Because of the large dipole of the free molecule (27 D in TB), electrostatic interactions play 
an important role in the LB film. We plot in Fig. [2] the intensity of the dipole, when we vary 
the molecular area in a free monolayer without a metallic electrode. At a molecular area of 
58.5 A 2 , the dipole is screened by a factor 2.6, due to the Coulomb interaction between the 
parallel dipoles. When the metallic electrodes are included, the dipole moment is enhanced, 
because of the attractive interaction between the charges in the dipole and those induced 
in the electrodes. We obtain a final dipole of 19 D in the diode at zero bias. This dipole 
layer gives rise to a built-in potential between the two sides of the film, as shown in Fig. [3], 
where we plot the average electrostatic potential in the system. In presence of the electrodes 
and at zero bias, opposite charges appear on the electrodes (5.8xl0~ 03 nelectrons/ A 2 ), so 
that the induced potential drop exactly cancels the one created by the dipole layer. When 
a bias is applied to the diode, the HOMO and the LUMO energies become linear functions 
of the voltage (V), with the same slope (77) for the two levels (Fig. HI), which is another 
consequence of their derealization. Indeed, if the two states wre localized at opposite sides 
of the molecule, the energy of the state closer to the grounded electrode would vary less 
than the other. For example, it was shown recently on D-cr-A molecules that the HOMO- 
LUMO gap is strongly distorted by the electric fielcr^, because the two states are localized 
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on different sites. In contrast, in our case of a D-7T-A molecule, the gap is almost independent 
of the applied voltage. The slope rj, obtained from the self-consistent calculation, strongly 
depends on the length of the alkyl chain (Fig. E]). It can be estimated, by comparing the 
center of the chromophore of the molecule with the molecular length. Thus, for Q-3CNQ, 
r] = 0.49 (close to 1/2), because the potential drop is symmetrically shared in the molecule, 
while for molecule 1 it is 0.21 (close to 1/4), because the states are localized in the Q-3CNQ 
part, close to the grounded electrode, whereas a large part of the potential drop takes place 
in the alkyl chains, which have a small polarizability due to their large gap^. This important 
effect is particularly striking in Fig. [31 where we compare the electrostatic potential in the 
system at zero bias and at +2 V. The voltage drop in the molecular layer, and thus the 
factor rj, can be simulated, using a simple dielectric model, where the different parts of the 
structure are described by a dielectric medium. We consider three regions: the vacuum gaps 
between the molecules and the metal electrodes, the Q-3CNQ part, and the alkyl chain, 
with respective widths d vacuumi d^, d a ik y i, and respective dielectric constants e n and e a ik y i- 
We obtain easily: 

2 j &7T d a lkyl 

alkyl- ((In -\- dvaccuum ' ^7r) 

This equation explains well the dependence of rj with the length of the alkyl chain (Fig. 
[5]). The best fit is obtained with = 10.0 and e a ik y i = 1.6, in agreement with the fact that 
the alkyl chains have a small polarizability. Measurements on monolayers of alkyl chains 
between metal electrodes give values of e a ik y i between 2.0 and 2.5, but for a density of one 
molecule per surface of 20-22 A 2 — (compared to 58.5 A 2 in our case). Our calculations for 
such a monolayer with a molecular area of 20 A 2 gives e a ik y i = 2.1 , in agreement with the 
experiments. 

III. CALCULATION OF THE CURRENT 

The current intensity is calculated using the Landauer formula [18] 

1= { J Lfz(e) - f R (e)}Tr{G%e)T R (e)G r (e)r L (e)}de (3) 
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where and fnie) are the Fermi-Dirac distribution functions in the left and right 

leads, respectively. G r ' a (e) are the retarded and advanced Green's functions of the molecule, 



and the T L,R (e) matrices that describe the coupling of the molecule to the electrodes 4 ". 
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The Green's functions of the semi-infinite metal electrodes are calculated in TBi*, using 
the decimation method 2 ^. The hopping matrix elements between metal orbitals (7) and 
molecular orbitals (f3) are restricted to the interactions with metal atoms of the top surfaces. 
They are written 

where d is the inter-atomic distance, do is the sum of the covalent radii of the two atoms, 
and V3(do) is the hopping integral, calculated from Harrison's rules^- at the inter-atomic 
distance is the hopping matrix element, as calculated by Client to study the current in a 
scanning tunneling microscope. The advantages of expression are: i) at small distance 
d close to do, the couplings are close to Harrison's values, which give a good description of 
the chemical interactions 21 ; ii) at larger distances (d >> do), the exponential dependence of 
the hopping integrals with d is correctly described by the expressions given by Chen^. For 
a given bias, we calculate the self-consistent electronic structure and the Green's functions 
of the molecule in the LB film. The Green's functions of the whole system are calculated by 
coupling the molecules to the leads. 



IV. RESULTS FOR AL ELECTRODES 



Let us consider the case of Al electrodes. We need to discuss the problem of the position 
of the molecular levels with respect to the Fermi level of the electrodes at zero bias. We have 
considered three situations (a, b, c), which are summarized in Fig. El The first situation 
(a) corresponds to the vacuum level alignment, a common approximation for physisorbed 
molecules. The HOMO and LUMO energy levels are fixed, with respect to the metal work 
function (4.2 eV), by the ionization potential I and the electron affinity A, respectively. 
A and I include the charging energy of the molecule, corresponding to its ionization when 
the current flows through the structure. We have calculated these values in LDA (J = 6.9 
eV, A = 2.7 eV for 1 and 2). But, for an ionized molecule in the diode, we must also 
consider the screening due to the electrodes, i.e. the interaction energy between the extra 
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charge (an electron in the LUMO if) e or a hole in the HOMO iph and the image charges 
induced in the electrodes 24 . It is given fairly accurately, in first-order perturbation theory, 
by E/i( e ) =< iih(e)\Vind\iJh(e) >/2 , where Vmd is the electrostatic energy of a charge between 
two metallic electrodes. We obtain in TB S e = 0.53ey and = 0.62eV for 1 (resp. 0.90 
eV and 0.97 eV for 2). We fix the molecular levels by applying a rigid shift to all unoccupied 
(occupied) states, such that the LUMO (HOMO) level position is given by A + E e (I- ). 
Fig. [6] (case a) shows that the Al Fermi level is closer to the LUMO level than the HOMO. 
In fact, the vacuum level alignment does not occur in many metal-organic systems^. In our 
case, the alignment may be influenced by charge distributions in the thin Al oxide, which 
is unavoidably present between the Al electrodes and the LB monolayer-. The main effect 
is a rigid shift of the whole set of molecular energy levels^. Thus, in order to investigate a 
broad range of possible experimental situations, we have considered a second situation (b) 
corresponding to the Al Fermi level at the middle of the HOMO-LUMO gap and, a third 
situation (c), symmetric to the first one, where the Al Fermi level is closer to the LUMO 
level. The calculated I(V) characteristics are plotted in Fig. [7] for molecule 1, with its long 
alkyl chain. In the level configuration (a), we predict an asymmetric /(V) curve with a 
forward behavior at negative bias, contrary to the Aviram-Ratner situation. The onset of 
the current at -1.1 V corresponds to the resonance of the LUMO level with the Fermi level 
of the D(+)-side electrode, as shown in Fig. HI The electrical rectification arises from the 
asymmetric profile of the electrostatic potential across the system, leading to rj ^ ~. The 
resonance of the LUMO with the Fermi level of the A(-)-side electrode could only occur at 
high positive bias (3.7 V). In configuration (b), the resonance with the LUMO shifts to a 
more negative voltage (-1.8 V) but there is an additional onset at +1.8 V, corresponding 
to the resonance with the HOMO level (almost invisible in Fig. [7] because the transmission 
coefficient is small). Finally, in configuration (c), only this latter resonance is possible at +1.0 
V, giving rise to a rectification at positive bias. The situation is completely different in the 
case of Q-3CNQ (inset of Fig. |7J), where rj ~ ~. If a Fermi level of an electrode energetically 
reaches some molecular level at a given bias, the same effect occurs symmetrically with the 
Fermi level of the other electrode at the opposite bias (Fig. H|). This leads to -almost 26 - 
symmetric I(V) curves for the three configurations (a, b, c). Experimental studies^- have 
been made on molecules 1 with a long alkyl chain; some Al| monolayer |A1 devices show (i) 
a rectification at positive bias and relatively low current^, (ii) some have symmetrical I{V) 
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curves, and (iii) some have relatively high currents and rectification at negative bias^. The 
calculations presented here would make case (i) correspond to configuration c, where the 
HOMO level is close to the metal Fermi level, and case (iii) correspond to configuration 
a, where the LUMO level crosses the Fermi level of the Al electrode with bias applied to 
it. All this experimental variability is either due to molecules turning upside down after 
film deposition and before measurement^, or it is due to non-uniform charges in the thin Al 
oxide 5 . 

V. RESULTS FOR AU ELECTRODES 

Recent experiments^ have been performed on LB films of the same molecules but using 
oxide-free Au electrodes. All the operating devices (not short-circuited) are rectifying at 
positive bias with a threshold at ~ +1.6 V (Fig. [S]). The main difference between Au 
and Al electrodes comes from the metal work function (resp. 5.3 eV and 4.2 eV). Shifting 
the zero bias Fermi level at -5.3 eV, we can just transpose the results of Fig. HJ from 
which we calculate the voltages for the resonance between the HOMO-LUMO levels and 
the Fermi levels. As oxide charges cannot be invoked in the case of Au electrodes, we 
have only considered the situation (a) where levels are aligned with respect to the vacuum 
level. We obtain in this case a rectification at positive bias (+1.24 V) in agreement with 
the experiments, corresponding to the resonance of the HOMO with the Fermi level at the 
D(+)-side electrode (situation close to the case c for Al). The resonance with the Fermi level 
at the A(-)-side only occurs at -4.6 V. We compare in Fig. [8] one typical experimental I(V) 
curved with the calculated one for Au|Ci6H 3 3Q-3CNQ|Au. As discussed above, the theory 
predicts a rectification at positive bias, in agreement with the experiments. But we have 
been obliged to multiply the calculated current by a factor 10 3 to be on the experimental 
scale. We have checked that there is no way to explain this factor by reducing the distance 
between the molecules and the electrodes. Thus one possible explanation is that some metal 
of the top Au electrode deposited on the LB film has slightly diffused in the molecular 
layer. This is equivalent to a reduction of the length of the alkyl chain, which acts as a 
barrier for the tunneling of the electrons. Fig. [8] shows that a correct magnitude of the 
current is obtained for Au|C 8 H 17 Q-3CNQ|Au, or equivalently for a diffusion of gold over a 
length of ~ 4 A(taking into account the tilting of the molecules). This is quite realistic with 



gold electrodes, even if the experimental procedure of Ref.-S has been designed to avoid at 
maximum the diffusion of the metal. The theory predicts a saturation of the current at high 
voltage, which is not observed experimentally as the experimental threshold (~ +1.6 V) 
is at higher voltage. We must add that the step-like dependence of the calculated current 
would be broadened by taking into account the interactions between the molecules or the 
coupling to molecular vibrations. On the experimental side, the application of high bias 
starts to degrade the molecular diode^. 

VI. CONCLUSION 

In conclusion, we have presented detailed calculations of the I(V) characteristics in a LB 
film of C16H33Q-3CNQ sandwiched between two Al or Au electrodes. We show the impor- 
tance to include Coulomb interactions and screening effects in a self-consistent manner. In- 
escapable conclusions are that the placement of the electroactive part of the molecule within 
the gap between metal electrodes is very important, and that an important experimental 
issue in the future will be to control the band alignment at the organic-metal interface. 
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FIG. 1: Top (a): Schematic representation of C16H33Q-3CNQ and of the diode. Bottom (b): the 
Aviram-Ratner mechanism for molecular rectification. 

FIG. 2: Evolution of the dipole moment of C16H33Q-3CNQ in a monolayer without electrode as a 
function of the molecular area. 

FIG. 3: Electrostatic potential in the metal | C16H33Q-3CNQ film | metal system at zero bias 
(dashed line) and at +2 V (straight line). The potential is defined as an average value in a lattice 
unit cell. 

FIG. 4: Evolution of the LUMO and HOMO energies (relative to the vacuum level) as a function 
of the applied bias (straight lines: C16H33Q-3CNQ, dotted lines: Q-3CNQ). At zero bias, the 
molecular levels and the metal Fermi levels are as shown for case a of Fig. [6J The crossings of the 
LUMO and HOMO energies with the Fermi levels (dashed lines) correspond to current thresholds 
in the I(V) characteristics. 

FIG. 5: Variation of r] with the length of the molecules C^^x+iQ-SCNQ (* full self-consistent tight 
binding calculation; straight line: simple dielectric model), r/ is the slope of the linear dependence 
of the HOMO and the LUMO energies with respect to the applied voltage. 

FIG. 6: Three possible schemes for the position of the molecular energy levels with respect to 
the Fermi level of the leads at zero bias. I and A are, respectively, the ionization energy and the 
electron affinity of the isolated molecules. and S e account for the reduction of the charging 
energy in the diode due to the screening by the electrodes. 

FIG. 7: Calculated I(V) characteristics of AIIC16H33Q-3CNQJAI for the three configurations of 
the molecular levels, with respect to the Fermi level of the electrodes, as defined in Fig. 6 (straight 
line: a, dashed line: b, dotted line: c). Insets: same for A1|Q-3CNQ|A1. 

FIGURES 
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FIG. 8: Experimental 6 and calculated current-voltage curves of molecular diodes with Au elec- 
trodes. 
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Abstract 

The current-voltage characteristics in Langmuir-Blodgett monolayers of 7-hexadecylquinolinium 
tricyanoquinodimethanide (C16H33Q-3CNQ) sandwiched between Al or Au electrodes is calculated, 
combining ab initio and self-consistent tight binding techniques. The rectification current depends 
on the position of the LUMO and HOMO relative to the Fermi levels of the electrodes as in the 
Aviram-Ratner mechanism, but also on the profile of the electrostatic potential which is extremely 
sensitive to where the electroactive part of the molecule lies in the monolayer. This second effect 
can produce rectification in the direction opposite to the Aviram-Ratner prediction. 
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I. INTRODUCTION 



Molecular electronics has gained large-scale interest in recent years. A promising sys- 
tem is the electrical rectifier based on suitably engineered molecules. At the origin of this 
idea, Aviram and Ratner 1 proposed in 1974 to use D-a-A molecules, where D and A are, 
respectively, an electron donor and an electron acceptor, and a is a covalent "sigma" bridge. 
Electrical rectification was observed only recently, in particular in Langmuir-Blodgett (LB) 
multilayers or monolayers of 7-hexadecylquinolinium tricyanoquinodimethanide (Ci 6 H 33 Q- 
3CNQ, 1, Fig. 1) sandwiched between metallic electrodes 2-6 . Even if these results represent 
an important progress to achieve molecular electronics, the physical mechanism responsible 
for the rectification is not clear. One critical issue is to know if the Aviram-Ratner model can 
be applied to 1, because it is a D(+)-7r-A(-) molecule 5 . On the theoretical side, these molec- 
ular diodes are complex systems, characterized by large and inhomogeneous electric fields, 
which result from the molecular dipoles 5 in the layer, the applied bias, and the screening 
induced by the molecules themselves and the metallic electrodes. A theoretical treatment of 
these effects is presently lacking, as it requires a self-consistent resolution of the quantum- 
mechanical problem, including the effect of the applied bias on the electronic structure. Our 
aim in this paper is to present such a theory, and its application to the systems experimen- 
tally studied in Refs. 3 ' 4 ' 6 . Combining ab initio and semi-empirical calculations, we show that 
the direction of easy current flow (rectification current) depends not only on the position 
of the HOMO and LUMO, relative to the Fermi levels of the metal electrodes before bias 
is applied, but also on their shift after the bias is applied: this situation is more complex 
than the Aviram-Ratner mechanism, and can provide a rectification current in the opposite 
direction. We calculate that the electrical rectification results from the asymmetric profile 
of the electrostatic potential across the system. We obtain that the Ci6H 33 tail plays an 
important role in this asymmetry, and we predict a more symmetric I(V) curve in the case 
of molecules with a small alkyl chain. Quite generally, our work emphasizes the importance 
of the electrostatic potential profile in a molecular system and suggests that this profile 
could be chemically designed to build new devices. 
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II. ELECTRONIC STRUCTURE AND POTENTIAL 



To calculate the electronic structure and the geometry of the free molecules (gas phase), 
we have performed ab initio calculations, using the DMOL code 7 in the local density approx- 
imation (LDA) with the functional of Ref. 8 , and in the generalized gradient approximation 
(GGA) of Ref. 9 . The two approximations give very close results. We use a double numeri- 
cal basis set (two atomic orbitals for each occupied orbital in the free atom), together with 
polarization functions (2p for H, and 3d for N and C). As the molecular diode is a complex 
system, which cannot be treated in LDA, we have developed self-consistent tight binding 
(TB) calculations, as described in Ref. 10 . We start from the GGA optimized geometries. C 
and N atoms are represented by one s and three p atomic orbitals, and H atoms by one s 
orbital. The diagonal terms of the Hamiltonian matrix are charge-dependent 10 ' 11 (in atomic 
units): 

H ia , ia = Hl ia - £ Qj - eV ext (i) (1) 

where a is the orbital index, i and j are atomic indices (j runs over all the atoms of 
the system), Qj is the net charge on atom j, Rij is the distance between atoms i and j, 
Uq is the intra-atomic Coulomb energy and the quantities Hf a ia define the s and p atomic 
levels. V ext (i) is the potential induced by the charges on the electrodes. It is calculated 
using a self-consistent solution of the Schrodinger and Poisson equations, with the boundary 
condition imposed by the applied bias. The second term in eqn. 1 is evaluated using the 
method of Ewald 12 , taking into account the two-dimensional periodicity of the LB film. The 
molecules are ordered on a 9Ax 6.5A rectangular lattice, corresponding to a molecular area 
of 58.5 A 2 , close to the experimental value of ~ 50 A 2 3 . The molecules are inclined by 45°, 
as measured experimentally 3 . Chemical interactions between molecules are supposed to be 
negligible. The LB film is sandwiched between two planar Al or Au electrodes, which we put 
at 3.2 A from the edges of the molecule, a reasonable value for molecules in weak interaction 
with a metallic surface 13 . We applied our calculations to 1, but 3 also to molecules with a 
smaller alkyl chain or without alkyl chain (Q-3CNQ, 2) for reasons which will become clear 
later. As expected, the electronic structure cose to the gap does not depend on the length 
of the alkyl chain for the free molecules. We predict a large dipole moment of 21 D in LDA 
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and 27 D in TB. The general results in LDA and in TB are quite similar, and agree with 
recent calculations 15 . Fig. 1 describes the original Aviram-Ratner mechanism 1 for a D-cr-A 
molecule. The electronic states are supposed to be totally localized either on the D side or 
on the A side. The HOMO(D) and LUMO(D) are high in energy, compared to, respectively, 
the HOMO (A) and LUMO(A). Therefore, a current can be established at relatively small 
positive bias, such that the Fermi level at the A side is higher than LUMO(A), and the 
Fermi level at the D side is lower than HOMO(D), providing that the electrons can tunnel 
inelastically through the a bridge. Thus we expect an asymmetric I(V) curve with an 
important rise of the current at positive bias above a given threshold (the electrode at A 
side is grounded, as in the experiments 3 ' 4 ' 6 ). In the case of our D(+)-7r-A(-) molecules 
(1 and 2), LDA and TB calculations show that the HOMO and the LUMO are strongly 
delocalized on the whole molecule. Even if the HOMO is more localized on the A part and 
the LUMO on the D part, the extension of the wave functions on the opposite side exceeds 
30%. Thus the two parts of the molecules cannot be decoupled, as in the original Aviram- 
Ratner model. This is obviously due to the n bridge between the donor and acceptor sites. 
Because of the large dipole of the free molecule (27 D in TB), electrostatic interactions play 
an important role in the LB film. We plot in Fig. 2 the intensity of the dipole, when we vary 
the molecular area in a free monolayer without a metallic electrode. At a molecular area of 
58.5 A 2 , the dipole is screened by a factor 2.6, due to the Coulomb interaction between the 
parallel dipoles. When the metallic electrodes are included, the dipole moment is enhanced, 
because of the attractive interaction between the charges in the dipole and those induced 
in the electrodes. We obtain a final dipole of 19 D in the diode at zero bias. This dipole 
layer gives rise to a built-in potential between the two sides of the film, as shown in Fig. 3, 
where we plot the average electrostatic potential in the system. In presence of the electrodes 
and at zero bias, opposite charges appear on the electrodes (5.8xl0~ 03 nelectrons/ A 2 ), so 
that the induced potential drop exactly cancels the one created by the dipole layer. When 
a bias is applied to the diode, the HOMO and the LUMO energies become linear functions 
of the voltage (V), with the same slope (rj) for the two levels (Fig. 4), which is another 
consequence of their derealization. Indeed, if the two states wre localized at opposite sides 
of the molecule, the energy of the state closer to the grounded electrode would vary less 
than the other. For example, it was shown recently on D-cr-A molecules that the HOMO- 
LUMO gap is strongly distorted by the electric field 15 , because the two states are localized 
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on different sites. In contrast, in our case of a D-7T-A molecule, the gap is almost independent 
of the applied voltage. The slope 77, obtained from the self-consistent calculation, strongly 
depends on the length of the alkyl chain (Fig. 5). It can be estimated, by comparing the 
center of the chromophore of the molecule with the molecular length. Thus, for Q-3CNQ, 
77 = 0.49 (close to 1/2), because the potential drop is symmetrically shared in the molecule, 
while for molecule 1 it is 0.21 (close to 1/4), because the states are localized in the Q-3CNQ 
part, close to the grounded electrode, whereas a large part of the potential drop takes place 
in the alkyl chains, which have a small polarizability due to their large gap 16 . This important 
effect is particularly striking in Fig. 3, where we compare the electrostatic potential in the 
system at zero bias and at +2 V. The voltage drop in the molecular layer, and thus the 
factor i], can be simulated, using a simple dielectric model, where the different parts of the 
structure are described by a dielectric medium. We consider three regions: the vacuum gaps 
between the molecules and the metal electrodes, the Q-3CNQ part, and the alkyl chain, 
with respective widths d vacuum , d,r, & a ik y u an d respective dielectric constants e n and e a ik y i- 
We obtain easily: 

77 = 1 ^ 7 — (2) 

2 U-alkyl 

alkyl- (d-n^ d vaccuurn • 6-^) 

This equation explains well the dependence of rj with the length of the alkyl chain (Fig. 
5). The best fit is obtained with e n = 10.0 and e a ik y i = 1-6, in agreement with the fact that 
the alkyl chains have a small polarizability. Measurements on monolayers of alkyl chains 
between metal electrodes give values of e a ik y i between 2.0 and 2.5, but for a density of one 
molecule per surface of 20-22 A 2 17 (compared to 58.5 A 2 in our case). Our calculations for 
such a monolayer with a molecular area of 20 A 2 gives e a ik y i = 2.1 , in agreement with the 
experiments. 

III. CALCULATION OF THE CURRENT 

The current intensity is calculated using the Landauer formula [18] 

1 = T I-J h{e) ~ ^(^M^r^G^r^H ( 3 ) 



5 



where /z,(e) and /ij(e) are the Fermi-Dirac distribution functions in the left and right 
leads, respectively. G r,a (e) are the retarded and advanced Green's functions of the molecule, 
and the T L ' R (e) matrices that describe the coupling of the molecule to the electrodes 18 . 
The Green's functions of the semi-infinite metal electrodes are calculated in TB 19 , using 
the decimation method 20 . The hopping matrix elements between metal orbitals (7) and 
molecular orbitals (/?) are restricted to the interactions with metal atoms of the top surfaces. 
They are written 

v* = |#^(4>) (4) 

where d is the inter-atomic distance, do is the sum of the covalent radii of the two atoms, 
and Vjp(do) is the hopping integral, calculated from Harrison's rules 21 at the inter-atomic 
distance is the hopping matrix element, as calculated by Chen 22 to study the current in a 
scanning tunneling microscope. The advantages of expression (4) are: i) at small distance 
d close to do, the couplings are close to Harrison's values, which give a good description of 
the chemical interactions 21 ; ii) at larger distances (d >> do), the exponential dependence of 
the hopping integrals with d is correctly described by the expressions given by Chen 22 . For 
a given bias, we calculate the self-consistent electronic structure and the Green's functions 
of the molecule in the LB film. The Green's functions of the whole system are calculated by 
coupling the molecules to the leads. 



IV. RESULTS FOR AL ELECTRODES 



Let us consider the case of Al electrodes. We need to discuss the problem of the position 
of the molecular levels with respect to the Fermi level of the electrodes at zero bias. We have 
considered three situations (a, b, c), which are summarized in Fig. 6. The first situation 
(a) corresponds to the vacuum level alignment, a common approximation for physisorbed 
molecules. The HOMO and LUMO energy levels are fixed, with respect to the metal work 
function (4.2 eV), by the ionization potential / and the electron affinity A, respectively. 
A and I include the charging energy of the molecule, corresponding to its ionization when 
the current flows through the structure. We have calculated these values in LDA (I = 6.9 
eV, A = 2.7 eV for 1 and 2). But, for an ionized molecule in the diode, we must also 
consider the screening due to the electrodes, i.e. the interaction energy between the extra 
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charge (an electron in the LUMO or a hole in the HOMO ip h and the image charges 
induced in the electrodes 24 . It is given fairly accurately, in first-order perturbation theory, 
by Sft(e) =< ifJh(e)\Vi n d\ifJh(e) >/2 , where Vi n d is the electrostatic energy of a charge between 
two metallic electrodes. We obtain in TB S e = 0.53eU and = 0.62eV for 1 (resp. 0.90 
eV and 0.97 eV for 2). We fix the molecular levels by applying a rigid shift to all unoccupied 
(occupied) states, such that the LUMO (HOMO) level position is given by A + S e (I- S ft ). 
Fig. 6 (case a) shows that the Al Fermi level is closer to the LUMO level than the HOMO. 
In fact, the vacuum level alignment does not occur in many metal-organic systems 25 . In our 
case, the alignment may be influenced by charge distributions in the thin Al oxide, which 
is unavoidably present between the Al electrodes and the LB monolayer 4 . The main effect 
is a rigid shift of the whole set of molecular energy levels 25 . Thus, in order to investigate a 
broad range of possible experimental situations, we have considered a second situation (b) 
corresponding to the Al Fermi level at the middle of the HOMO-LUMO gap and, a third 
situation (c), symmetric to the first one, where the Al Fermi level is closer to the LUMO 
level. The calculated I(V) characteristics are plotted in Fig. 7 for molecule 1, with its long 
alkyl chain. In the level configuration (a), we predict an asymmetric I(V) curve with a 
forward behavior at negative bias, contrary to the Aviram-Ratner situation. The onset of 
the current at -1.1 V corresponds to the resonance of the LUMO level with the Fermi level 
of the D(+)-side electrode, as shown in Fig. 4. The electrical rectification arises from the 
asymmetric profile of the electrostatic potential across the system, leading to rj ^ \. The 
resonance of the LUMO with the Fermi level of the A(-)-side electrode could only occur at 
high positive bias (3.7 V). In configuration (b), the resonance with the LUMO shifts to a 
more negative voltage (-1.8 V) but there is an additional onset at +1.8 V, corresponding 
to the resonance with the HOMO level (almost invisible in Fig. 7 because the transmission 
coefficient is small). Finally, in configuration (c), only this latter resonance is possible at +1.0 
V, giving rise to a rectification at positive bias. The situation is completely different in the 
case of Q-3CNQ (inset of Fig. 7), where r\ « \. If a Fermi level of an electrode energetically 
reaches some molecular level at a given bias, the same effect occurs symmetrically with the 
Fermi level of the other electrode at the opposite bias (Fig. 4). This leads to -almost 26 - 
symmetric I(V) curves for the three configurations (a, b, c). Experimental studies 3-5 have 
been made on molecules 1 with a long alkyl chain; some Al| monolayer |A1 devices show (i) 
a rectification at positive bias and relatively low current 5 , (ii) some have symmetrical /(V) 
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curves, and (iii) some have relatively high currents and rectification at negative bias 4 . The 
calculations presented here would make case (i) correspond to configuration c, where the 
HOMO level is close to the metal Fermi level, and case (iii) correspond to configuration 
a, where the LUMO level crosses the Fermi level of the Al electrode with bias applied to 
it. All this experimental variability is either due to molecules turning upside down after 
film deposition and before measurement 3 , or it is due to non- uniform charges in the thin Al 
oxide 5 . 

V. RESULTS FOR AU ELECTRODES 

Recent experiments 6 have been performed on LB films of the same molecules but using 
oxide-free Au electrodes. All the operating devices (not short-circuited) are rectifying at 
positive bias with a threshold at ~ +1.6 V (Fig. 8). The main difference between Au 
and Al electrodes comes from the metal work function (resp. 5.3 eV and 4.2 eV). Shifting 
the zero bias Fermi level at -5.3 eV, we can just transpose the results of Fig. 4, from 
which we calculate the voltages for the resonance between the HOMO-LUMO levels and 
the Fermi levels. As oxide charges cannot be invoked in the case of Au electrodes, we 
have only considered the situation (a) where levels are aligned with respect to the vacuum 
level. We obtain in this case a rectification at positive bias (+1.24 V) in agreement with 
the experiments, corresponding to the resonance of the HOMO with the Fermi level at the 
D(+)-side electrode (situation close to the case c for Al). The resonance with the Fermi level 
at the A(-)-side only occurs at -4.6 V. We compare in Fig. 8 one typical experimental I(V) 
curve 6 with the calculated one for Au|Ci 6 H 33 Q-3CNQ|Au. As discussed above, the theory 
predicts a rectification at positive bias, in agreement with the experiments. But we have 
been obliged to multiply the calculated current by a factor 10 3 to be on the experimental 
scale. We have checked that there is no way to explain this factor by reducing the distance 
between the molecules and the electrodes. Thus one possible explanation is that some metal 
of the top Au electrode deposited on the LB film has slightly diffused in the molecular 
layer. This is equivalent to a reduction of the length of the alkyl chain, which acts as a 
barrier for the tunneling of the electrons. Fig. 8 shows that a correct magnitude of the 
current is obtained for Au|CsHi7Q-3CNQ|Au, or equivalently for a diffusion of gold over a 
length of ~ 4 A(taking into account the tilting of the molecules). This is quite realistic with 
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gold electrodes, even if the experimental procedure of Ref. 6 has been designed to avoid at 
maximum the diffusion of the metal. The theory predicts a saturation of the current at high 
voltage, which is not observed experimentally as the experimental threshold (~ +1.6 V) 
is at higher voltage. We must add that the step-like dependence of the calculated current 
would be broadened by taking into account the interactions between the molecules or the 
coupling to molecular vibrations. On the experimental side, the application of high bias 
starts to degrade the molecular diode 6 . 

VI. CONCLUSION 

In conclusion, we have presented detailed calculations of the I(V) characteristics in a LB 
film of Ci6H 33 Q-3CNQ sandwiched between two Al or Au electrodes. We show the impor- 
tance to include Coulomb interactions and screening effects in a self-consistent manner. In- 
escapable conclusions are that the placement of the electroactive part of the molecule within 
the gap between metal electrodes is very important, and that an important experimental 
issue in the future will be to control the band alignment at the organic-metal interface. 
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FIG. 1: Top (a): Schematic representation of C16H33Q-3CNQ and of the diode. Bottom (b): the 
Aviram-Ratner mechanism for molecular rectification. 

FIG. 2: Evolution of the dipole moment of C16H33Q-3CNQ in a monolayer without electrode as a 
function of the molecular area. 

FIG. 3: Electrostatic potential in the metal | C16H33Q-3CNQ film | metal system at zero bias 
(dashed line) and at +2 V (straight line). The potential is defined as an average value in a lattice 
unit cell. 

FIG. 4: Evolution of the LUMO and HOMO energies (relative to the vacuum level) as a function 
of the applied bias (straight lines: C16H33Q-3CNQ, dotted lines: Q-3CNQ). At zero bias, the 
molecular levels and the metal Fermi levels are as shown for case a of Fig. 6. The crossings of the 
LUMO and HOMO energies with the Fermi levels (dashed lines) correspond to current thresholds 
in the I(V) characteristics. 

FIG. 5: Variation of rj with the length of the molecules C x H2 X +iQ-3CNQ (* full self-consistent tight 
binding calculation; straight line: simple dielectric model), rj is the slope of the linear dependence 
of the HOMO and the LUMO energies with respect to the applied voltage. 

FIG. 6: Three possible schemes for the position of the molecular energy levels with respect to 
the Fermi level of the leads at zero bias. I and A are, respectively, the ionization energy and the 
electron affinity of the isolated molecules. and S e account for the reduction of the charging 
energy in the diode due to the screening by the electrodes. 

FIG. 7: Calculated I(V) characteristics of AIIC16H33Q-3CNQJAI for the three configurations of the 
molecular levels, with respect to the Fermi level of the electrodes, as defined in Fig. 6 (straight 
line: a, dashed line: b, dotted line: c). Insets: same for A1|Q-3CNQ|A1. 

FIGURES 
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FIG. 8: Experimental 6 and calculated current-voltage curves of molecular diodes with Au elec- 
trodes. 
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